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We study superpotential perturbations of q deformed = 4 Yang-Mills for q a 
root of unity. This is a special case whose geometry is associated to an orbifold 
with three lines of codimension two singularities meeting at the origin. We perform 
field theory perturbations that leave only co-dimension three singularities of conifold 
type in the geometry. We show that there are two "fractional brane" solutions of 
the F-term equations for each singularity in the deformed geometry, and that the 
number of complex deformations of that geometry also matches the number of sin- 
gularities. This proves that for this case there are no local or non-local obstructions 
to deformation. We also show that the associated Dijkgraaf-Vafa matrix model has 
a solvable sector, and that the loop equations in this sector encode the full deformed 
geometry of the theory. 
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I. INTRODUCTION 



Supersymmetric field theories are a very attractive area of modern theoretical physics 
not only because of their possible application to physics beyond the standard model, but 
also because they have features that can be solved exactly non-perturbatively Q]. The 
program of geometric engineering produces these type of field theories by placing compact 
D-brane configurations on non-compact Calabi-Yau manifolds. Qne of the main features of 
using non-compact CY manifolds is that four dimensional gravity decouples, and in the low 
energy limit we are left with an ordinary supersymmetric field theory. The advantage of 
such a formulation is that many aspects of the field theory can be understood geometrically. 
For example, the holomorphic structure of the Calabi-Yau manifold might be deformed at 
singularities, and this can be associated to confinement 0, 0, |^ 

Having this geometric information is important in terms of the AdS / CFT correspondence 
if we want to understand gravity duals to certain field theories, and it is also useful for model 
building if we want to have a theory on a CY throat that describes the standard model of 
particle physics (or some similar theory) via a D-brane configuration. In practice it is very 
hard to compute the low energy field theory on the D-brane because we need to have full 
knowledge of the superpotential, and this can involve calculating arbitrarily complicated n 
point functions on the D-brane. This can be overcome if we have a large amount of symmetry 
that forbids most terms (this is what happens in practice for CFT's associated to toric CY 



singularities, see for example Q ). 

We can also try to work in the inverse direction. Start with some field theory and try to 
construct a dual geometry that realizes it. Some progress in this direction has been made in 
y^Ql; by turning the problem into calculating certain properties and features of an algebra 
associated to the field theory. However, there is no algorithm to systematically carry out 
the procedure outlined there (except when there is a lot of symmetry), even at the classical 
level, and one has to work with specific examples. 

One of the big discoveries of the geometric engineering program was the realization by 
Dijkgraaf and Vafa that one can calculate an effective superpotential for gaugino conden- 
sation by studying the planar diagrams of an associated matrix model derived from the 
superpotential y, Q . This has been derived exactly in field theory in the works 11. 1^. 
In many examples, one only has to worry about a one-matrix model, and one can relate the 
spectral curve of the matrix model to the (deformed) geometry of the original Calabi-Yau 
space. This usually works best if one can understand the model as a deformation of an 
Af = 2 theory, where one also can solve the model using Seiberg-Witten techniques ^| . 

In most cases however, one has to try and understand the problem in terms of a more 
generic looking matrix model. In this case it is not obvious that there is a spectral curve, 
or how the geometry of the CY manifold is related to the matrix model. A proposal for 



what one is supposed to do in this case has been outlined in by studying the quantum 
moduli space of a brane in the bulk, but to date there is no proof that this always works. 
One can also encode in some examples the CY complex structure geometry by looking at 
the loop equations associated to the matrix model [3] , see also 13, ll^ for other examples. 

The purpose of this paper is to explore this issue for a particular class of models. The 
strategy that one can adopt is to begin in a configuration where one knows the field theory 
very precisely. Then one has a dictionary between certain superpotential deformations and 
certain deformations of the geometry. This is usually simple to understand in the case where 
one has (locally) an ADE orbifold singularity (these are codimension two singularities in the 
CY geometry) and one can identify how the closed string moduli that resolve the singularity 
couple to a D-brane. In this case, one can follow the deformation at the D-brane level (field 
theory) as well and one can try to reconcile both geometries. In our examples we begin with 
the C^/Z„ X Z„ orbifold in the presence of discrete torsion and follow that general idea in 
the example. A generic deformation will resolve the codimension two singularities and might 
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leave behind some co dimension three singularities in the geometry. Most of the time, these 
are just conifold singularities (these are all the singularities we find in the generic case), and 
we understand how these get deformed locally by placing fractional branes at the singularity 

We expect this behavior to be universal in the presence of fractional branes at the sin- 
gularity and to have one deformation parameter per conifold singularity. However, there 
always exist the possibility that some of these deformations are obstructed (perhaps non- 
locally), or that the naive geometric intuition does not match the field theory analysis. 
Indeed, examples have been found where there is an inconsistency of the chiral ring once 
gaugino condensates are taken into account [2^ |^, and one finds that supersymmetry is 



broken or there is runaway behavior [22|. Morever, even if we understand the local behavior 
at each singularity, we might have trouble understanding the global deformed geometry in 
detail. This is especially hard in the case of non-complete intersection varieties, although 
recent progress has been made for the special case of toric geometries j^]. In this paper we 
show that for the example we study there are no such obstructions and that the field theory 
analysis reproduces the geometric analysis precisely. 

The model we work with is built by starting with the q deformed J\f = 4 SYM theory, 
which is associated to the C^/Z„ x Z„ orbifold with discrete torsion 0, [2^. The theory 
has the same matter content of A/" = 4 SYM, but a different superpotential 



with g" = 1. T 



W = tr{XYZ - qXYZ) (1) 
lis model has been studied extensively in the past, in various contexts, see 



for example |2Sl, |22l, 12^1, |22i , and more recently the interest in the model for general q has 
been sparked by the dual supergravity solution found by Maldacena and Lunin . All the 
models for different roots of unity give rise in the Dijkgraaf-Vafa program to three-matrix 
models which on the face look very similar to each other. The associated geometries end up 
being fairly different, although one can still study them in a unified context. 

Our reasons for studying this particular model are various. The theory has a unique gauge 



^ These brane configurations can also produce duality cascades in the field theory for the low energy 
dynamics of branes at the singularity. It also seems to be the case that the field theory associated 
to conifolds is always in the same universality class ,183. Recently, the understanding of the full moduli 
space of vacua of these cascading theories has been improved substantially |l9l | 
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group with a simple matter content. The corresponding matrix model has three matrices 
and looks simple as well. On the other hand the geometric interpretation is relatively 
complicated and leads to a very useful laboratory to study complex structure deformation 
problems. 

The idea is now to deform the field theory superpotential given above by twisted sector 
deformations. These deformations are generated by 



for k not a multiple of n 1255. These are non-zero elements of the chiral ring whose support is 
concentrated at the codimension two singularities of the C^/Z„ x Z„ orbifold space. Because 
these are elements of the chiral ring with support on the singularities, these are associated 
to complex structure deformation moduli in string theory. This is what one obtains from 
an analysis within the framework of the AdS/CFT correspondence Isil Is^. 

However, we have to deal with the problem that the codimension two singularities are 
non-compact and there are an infinite number of possibilities, some of which might deform 
the geometry severely. This same problem is present in the work on deformations of = 2 
theories . The idea is to think of these deformations we are performing as being marginal 
or relevant. This depends on a notion of i?-charge, that might differ from the naive one 
that one deduces directly from the q deformed Af = 4 theory, as if one had a different UV 
fixed point. It must also be the case that since we want the theory to be dominated by the 
A/" = 4 g deformed superpotential, that this represents a marginal deformation. Obviously 
we have a U{1)^ symmetry of the q deformed potential. With this symmetry we can assign 
different R-charges to X, Y, Z so that they are all positive and given by 71, 72, 73 respectively, 
and so that the R symmetry of the q deformed superpotential is exactly two. This type of 
modification is similar to what is done in for non-anticommutative power counting. 
Here, it modifies the notion of relevant operators for superpotential couplings, which are 
the ones that control the holomorphic data of the theory. This in the end tuns out to be a 
convenient tool to keep track of holomorphy without doing a complicated analysis of charges. 

With this prescription, we are only allowed to add monomials tr(X^) such that k'ji < 2 
tr(r') so that /72 < 2 and tr(Z™) so that rri'js < 2. 

Indeed, we can choose 73 arbitrarily close to one, so we can deform both in the X and 
Y singularity as much as we want to, and then the only extra term that depends on Z is 



tr(X^),tr(r^),tr(Z'=) 



(2) 
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exactly tr(Z). For the sake of simplicity, we will take the coupling associated to tr(Z) to 
be exactly zero. In principle one can repeat the analysis we have done here for this slightly 
more general case. Part of this analysis at the classical level for a particular case has been 
done in [si] however, we have not yet found a closed form expression of the geometry for 
the full theory in the more general case. 

The general geometry that we get without gaugino condensations is described by the 
following simple hypersurface equation 

uvw = Q{t,u) +Q'{t,v) -r (3) 

where the orbifold with discrete torsion corresponds to the case uvw = t". Q and Q' are 
polynomials with leading term t", in fact it will turn out that the polynomials we find can 
be written in factorized form 

n-1 

Q{t,u) = l[{t-f{q^u'/n) (4) 

1=0 
n-1 

Q'{t,v) = l[it-g{qV/n) (5) 

i=0 

where on the right hand side we have to choose some n-th root of u,v and where /(O) = 
g{0) = 0. It is easy to see that it does not matter which root we choose, we get the same end 
result, and that Q and Q' are indeed polynomials of u, v. It is clear that Q has singularities 
when different roots of Q at m fixed coincide. The same is true for Q'. These give rise to 
singularities of the full geometry when either m or f are zero. These are generically conifold 
singularities. It is easy to convince oneself that these are the only singularities in the above 
geometry. The derivative with respect to w forces either m or w to be zero. After this is 
used, one is left with either Q or Q' only. The derivative with respect to the variable that 
is turned to zero is used in turn to determine the value of w. 

The main issue we will be concerned with is to count all the conifold points of this 
geometry. We will show that for each of them there is a log-normalizable deformation of the 
geometry in a suitable sense, and we will show that for each of these singularities one also 
finds a pair of fractional branes that can get trapped at the singularity. We furthermore 
show that if we take into account the loop equations of the associated matrix model, that we 
can recover the full geometry, and that there are no field theory obstructions to deformation. 
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This fits well with the conjectures in where it is stated that only in the case of 
geometric obstructed deformations we should get an inconsistency of the quantum deformed 
chiral ring of the theory. 

In a sense, our results are exactly as they should be and provide a technically involved 
consistency check of the Dijkgraaf-Vafa ideas, as well as an interesting solution to a com- 
plicated multi-matrix model. Presumably these techniques can be extended to other setups 
with multiple gauge groups as well. A very interesting case to analyze is a perturbation 
of the C^/Z2 X Z2 singularity without discrete torsion. In that case it has been observed 
that fractional branes at the origin trigger runaway behavior for fractional branes on the 
codimension two singularities ||^. Maybe in this case one would obtain a more involved 
obstruction to deformation. 

The paper is organized as follows: In section 2 we will analyze the deformed Af = 4 theory 
by solving the F-term constraints to find solutions corresponding to branes in the bulk as 
well as fractional brane solutions at singularities, and find the hypersurface equation defining 
the deformed geometry. We will also count the number of relevant quantum deformations, 
and show that it is equal to the number of singularities of the geometry. In section 3 we will 
analyze the theory by finding the quantum relations of the associated matrix model, and 
show how the loop equations encode the complex structure of the CY geometry. 



II. DEFORMED Af = A SUPERPOTENTIAL 

The unique A/" = 4 supersymmetric gauge theory in four dimensions has three adjoint 
fields X, Y, Z in the gauge group G and a superpotential 

ti{XYZ - XZY). (6) 

For our purposes G will be SU{N). We are interested in studying the structure of closely 
related field theories with the same field content, but with a different superpotential. Some 
of these are very interesting models in geometric engineering because they correspond to 
the field theories associated to D-branes in singular geometries. Some of these are orbifolds 
with discrete torsion C^/Z„ x Z„ 24,|27|. For this case the superpotential is: 



Wo = tT{XYZ - qXZY) 



(7) 
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where = 1 is a root of unity. The geometry is given by the hypersurface in defined by 

uvw = (8) 

and has three fines of singularities (along tfie X, Y, and Z axis) of codimension two wfiicfi 
intersect at tfie origin. At a generic point of tfie singular locus tfie geometry looks locally like 
a C^/Z„ X C quotient and tfie low energy effective field tfieory of a brane at tfiat locus fias 
an A/" = 2 supersymmetry. Globally tfiis is not true, but one can analyze tfie tfieory in tfie 
same spirit of Seiberg-Witten. Tfiis task fias been performed in a series of papers, mostly in 



tfie work of Dorey et al [36|. Tfie relation to AdS/CFT fias also been extensively studied, 
mostly because tfiese tfieories belong to a one parameter family of marginal deformations of 
A/" = 4 SYM [2^ and as sucfi one could expect tfiat for sufficiently small deformations one 
can understand tfie dual AdS space at tfie level of a supergravity analysis P, 21, 0,0- 

Perfiaps tfie biggest surprise in tfie study of tfiis class of deformations is tfiat tfie field 
tfieory is so ricfi, in spite of looking deceptively close to A/" = 4 SYM. We wisfi to furtfier de- 
form tfie geometry to break tfiese lines of singularities down to points, i.e. codimension tfiree 
singularities. Tfiis is accomplisfied by turning on additional superpotential term couplings. 

We want to argue tfiat in some sense tfiese are relevant deformations of tfie superpotential, 
and tfiat tfie superpotential Wo dominates tfie geometry. Tfiis can only be argued if tfie 
superpotential deformation satisfies some additional constraints; tfie ones tfiat are relevant 



by power counting fiave been analyzed in 



24 



27 



Tfiere is also tfie Polcfiinski-Strassler 



setup j38[, tfiat does not fall into tfie category of deformations tfiat we allow for g = 1. 

Here we will give a definition of wfiat relevant terms mean in a fiolomorpfiic sense, from 
tfie point of view of deformation tfieory. Tfie idea is simple, tfie fiypersurface equation 
uvw = f fias a U{lY symmetry under independent rescalings of u,v,w {t gets rescaled by 
eacfi of tfiese transformations). We can cfioose from tfiis set a generalized degree for m, f , w 
nri,nr2,nr3 sucfi tfiat tfiey all fiave positive degree, and sucfi tfiat tfie degree of t is one. 
In tfiis way ri + r2 + r-^ = 1. A complex structure deformation of tfie equation will be 
given by a polynomial / in tfie variables u,v,w,t, so tfiat tfie fiypersurface equation reads 
uvw = + f. Tfie deformation associated to / will be relevant if all tfie monomials of / 
fiave strictly smaller degree tfian n. At tfie level of tfie field tfieory, X, Y, Z will fiave cfiarge 
1^1,1^2, fs witfi respect to tfiis i?-like symmetry, and tfie deformation of tfie superpotential will 
be a polynomial in X, Y, Z of degree less tfian one. Tfie standard cfioice is ri = r2 = rs = 1/3. 
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We will consider a different choice where is arbitrarily close to 1. In this sense, we can 
allow any deformation of the superpotential which is a polynomial in X, Y, but which does 
not depend on Z. We will also not consider a linear term in Z because this simplifies the 
analysis. 



A. F-term equations 

We start from the most general form for a deformed superpotential with three fields 
X, y, Z in the adjoint: 

W = tr{XYZ - qXZY + V{X) + U{Y)) (9) 

where 

q'' = 1 (10) 
V{x)^ E J^' (11) 

i^O mod n 

u{y)- E jy' (12) 

i^O mod n 

The restrictions on V and U guarantee that this superpotential has a geometric brane 
realization. This is because the generic complex deformation we are turning on should 
correspond to twisted sector moduli localized at the codimension two singularities (these 
are the ones that resolve the singularities). The terms X^F™ have Z„ x Z„ charge {k,m) 
mod (n) under the quantum symmetry of the orbifold. Although the form above seems 
non-generic, note that any cross terms (with non-trivial twisted sector charge) of the form 
tr(X'^y"*) are proportional to a polynomial in X or y plus terms proportional to the F-term 
equations of motion for the field Z. These can be absorbed in a redefinition of Z by a linear 
shift. These conditions also make it so that the F-term equations are solvable in terms of 
n X n matrices (this corresponds to a brane in the bulk in the underformed case), and the 
moduli space of the brane in the bulk is exactly the deformed geometry. We will now derive 
this from studying the F term constraints. We have three matrix equations derived from 
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the F-term constraints: 



XY - qYX 
ZX - qXZ 

YZ - qZY 







i^O mod n 

E 

i^O mod n 



i-1 



i-1 



(13) 
(14) 

(15) 



and we want to find all solutions of these matrix equations, up to equivalence. This is the 
problem of solving the classical moduli space of vacua of the theory. 

In general, if we have two different solutions, we can take direct sums, and they also 
solve the equations. The most interesting solutions are those that are given by irreducible 



solutions (those that can not be written as direct sums o 
ones that we associate to individual branes following P,!?! 



smaller solutions. These are the 



B. Bulk Brane Solutions 

We first find the solutions to the F-term equations of dimensionality n, since those rep- 
resentations correspond to branes in the bulk. We begin by solving equation (fT^ . We can 
choose X to be diagonal, and then the general solution is of the form 



X = xP Y = yQ 



(16) 



where P and Q have the properties P" = = 1, PQ = qQP, and x and y are free complex 
parameters different from 0. When any of these two go to zero, one has to be more careful 
in describing the geometry, because either X" or Y" become zero. This means that all 
eigenvalues of X or y go to zero, but the matrices may still be upper or lower triangular, 
and this is necessary to be able to satisfy the equations (fT^ and (fT3j) 
With P chosen diagonal, the matrices are uniquely: 

/l \ /oo...oi\ 



T,n— 1 



Q 



1 ... 
1 ... 

\0 1 J 



(17) 
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the so called clock and shift matrices. Now we insert these into the equation ()14p : 

ZX-qXZ = - biv'-'Q'-^ (18) 

ij^O mod n 

It is a well known fact that the P and Q matrices span the n x n complex matrices as a 
vector space, so we can write Z in general a.s Z = 'Yll^=Q^ijQ^P^ ■ Inserting this into the 
equation allows us to solve it: 

n—X n—l 

CijxQ'P'^' - Y c^,xqPQ'P= = - b^y'^'Q''' (19) 

i,j=0 i,j=0 ij^O mod n 

n-1 

Y c.,xQ'P'^\l - q'^') = - Y b^y'-'Q'-' (20) 

i,j=0 iy^O mod n 

Considering the linear independence of the matrices, we can equate terms with equal powers 
of Q and P, giving for i + 1 7^ mod n: 

c,_i,„_ix(l - = - Y ^^y''' (21) 

j=i mod n 
j=i mod n 

For the other case, i + 1 = mod n: 

c„_i,,xP^(l-g") = (23) 
= (24) 

is satisfied for all c„_ij. Let us now relabel our sum and call Cj_i_n_i = {i runs from 1 to 
n — l). We can rewrite Z to make it a bit clearer: 

n— 1 n 

Z = Y ^Q'^'P""-' + Y Cn-i,Q"-'P^ (25) 

i=l i=0 

Now insert this form into equation p5|) 

YZ - qZY = - Y (26) 

27^0 mod n 

n—l n 

Y z[yQ'P''-\l - g") + J] c„_i,ygP^(l -q') = - Y (^^^'~'P"' (27) 

i=l i=0 i^O mod n 

n 

Y Cn-i,yP\l - q'^') = - Y <^i^'~'P''' (28) 

1=0 ijtO mod n 
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Exactly as in equation ()2ip above we equate terms with the same powers of P to get: 

Q 



cn-i,-i = - E (29) 

j=i mod n 



with i = l..n — 1, and that c„„i_„_i is unrestricted. We now make the definitions: c„_i^j_i = 
Zi and that c„_i^„_i = 2;. We can now write down the most general n dimensional solution 
to the F-term equations: 

X = xP (30) 

Y = yQ (31) 

n— 1 

Z = zQ'^-^P'^-^ + ^ {ziQ''-^P'~^ + z^Q'-^P"^-^) (32) 
where the Zj, z^' are given by 

E !'"'t^ (33) 



i=l 



q- 

j=i mod n 



j=j mod n 

This gives us a family of solutions parameterized by three complex numbers x, y, z, and 
therefore correspond to branes free to propagate in the bulk (which is a three dimensional 
CY geometry). However, the parameters x, y, z are not gauge invariant, as we have implicitly 
made gauge transformations by selecting X diagonal. One can show that x transforms as 
X ^ qx hj conjugating the solution with Q. This also transforms z and the Zi, z[, so the 
parameters x, y, z give a branched covering of the moduli space. 

We want to find the moduli space more explicitly as a hypersurface equation in C^, where 
we can make contact with the original, undeformed geometry. 

To find the proper equation, we need to find the gauge invariant observables of the 
configurations, and this can be translated into the center of the algebra defined by the three 
equations (the center is the set of all polynomials in X, F, Z which when evaluated in the 
solutions we found give matrices proportional to the identity.) This has been advocated as 
the correct description of D-branes at singularities in 

Note that X" = x" = det X = u and = ?/" = det Y = v are obviously in the center. 
It is also clear that is not in the center (which can be readily verified from the algebra). 
However, one can find that det Z = w is, by construction, in the center. (A note of caution: 
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w is not obviously a polynomial of X, Y, Z, but after some algebra we will find that it is 
). A final variable in the center is given by XYZ + f{X) + g{qY) = t, where / and g are 
polynomials. We can solve for them easily given the above solution: 

n-l 

XYZ ^xyz-Y, [xyz^P' + xyZiq'Q') (35) 

i=l 



i^O mod n i^O mod n 



SO 



/w= E r^^' (37) 



1 

ij^O mod n 

9iY)- E Y^l^^ (38) 

iy^O mod n 

Note that / and g are polynomials in X and y of the same order as V and U respectively. 
Also, we have that f{x) = xy^Zi and g{y) = xyY^ z\. 

An aside before we derive the hypersurface equation: How do we know that det Z is 
a polynomial in the algebra? We note that we can solve for the Zi^ in terms of f{q^x), 
specifically 

-^E^"W^)=^i (39) 
nxy ^ 

and likewise for the 2;^'s. We can use the characteristic equation for the Z matrix, which in 
general is of the form 

±A" + tr(Z) A"-^ + • • • + det Z = (40) 

where the omited terms are polynomial combinations of higher power traces. Note first 
that tr(Z) = 0, so there is no n — 1 power term. Plugging in the matrix Z itself into 
the characteristic equation (it must be a solution to its own characteristic equation by the 
Cayley-Hamilton theorem), we see that 

det Z = - polynomial in Z of degree n-1 (41) 

where the coefficients of the polynomial are combinations of the coefficients of Z we found 
above. As we showed, the coefficients of Z are simply the Zi% (and some g''s), which can be 
written in terms of polynomials of the X or y matrices (we must allow for dividing by x, y, 
but these are non-zero). It turns out that these terms in the denominators cancel. 
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We can continue on to obtain the full hypersurface equation. Using the expression for t, 
we simply take the determinant: 

XYZ = t - /(X) - g{qY) (42) 
uvw = <lei{t- f{X)~g{qY)) (43) 

This is relatively easy to analyze, since the t — f{X) is diagonal and the ~g{qY) is strictly 
off diagonal. This provides a factorization: 

n-1 

uvw = Y\_(^ " fi'f^)) + terms with g{y)'s (44) 

2=0 

We can go further. From the symmetry between X and Y present in the superpotential, we 
could have equally likely chosen Y = yP and X = xQ, so there must be symmetrical terms 
for each u and v. This gives: 

n— 1 n— 1 

UVW = Y\_(^ ~ fio"^)) + Wjy^ ^ gi'fy)) — t^ + cross terms with at least one / and g (45) 

From the original determinant expression, we also see that the RHS is invariant under 
X ^ qx and y qy {x qx implies f{X) Q^X^Q~^ for each term in f{X) which is 
invariant under the determinant, and y qy implies likewise with P instead). From this, we 
know that the right side can only be a function of the invariants under these transformations, 
namely u and v. Then from our further factoring, we can say that 

uvw = Q{t, u) + Q'{t, v)-f + uvR{t, u, v) (46) 

where R has degree less than n — 1 in all it's variables, and 

n-1 

Q{t,u) = l[{t-f{q'x)) (47) 

i=0 
n-1 

Q\t,v) = l[{t-g{q'y) (48) 

i=0 

We can now do a change of variables, shifting w ^ w + R{t,u,v), giving us our final 
hypersurface equation 

uvw = Q{t,u) +Q'{t,v) (49) 

which is surprisingly simple to write and analyze. When either g or f are zero it reduces to 
known geometries (see for example, j3|), and turning on both gives the simplest geometry 
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that one could possibly construct that captures both limits correctly. This is because any- 
thing more complicated that could appear on the right hand side could be shifted away by 
a coordinate redefinition of w. Notice that the only hard part of the analysis is to actually 
find w as a polynomial in X, F, Z, which we did implicitly. The most important thing is 
that in spite of beginning with some complicated looking solutions of the F-terms, the final 
form of the geometry is very simple. 



C. Fractional Brane Solutions 



The solutions we found above correspond to the generic bulk brane. It is interesting to 
find out if there are some cases where we can find solutions of smaller rank than n. In 
the case of orbifolds these solutions correspond to the existence of fractional branes at the 
singular locus. We will call these solutions the fractional brane solutions, and we will give 
them the standard interpretation, namely, as branes stuck at singularities of the geometry 
we found. Also, since it is easy to write down m, f , t in the geometry as polynomials of 
X, y, Z, and to use that information to calculate w, we will calculate these coordinates for 
the fractional brane solutions explicitly. 

Let us enumerate solutions with dimension p < n. If we begin by taking X to be diagonal, 
then equation (|T!^ gives: 



(50) 



with no sum implied. We want solutions where both X and Y are non-zero, which gives 
Uii = 0. We can also satisfy p — 1 of the rest of the equations by choices of the Xj's. We 
choose without loss of generality that qxi = Xj+i. Therefore only are non-zero, with 
the rest of the entries in Y zero. Relabeling xi = x and choosing a gauge where all the 
Ui-^i/s are 1, gives the following most general form for X and Y solving the first F-term 
equation ^ 



( 



X 



X 



1 ... 

q ... 








Y 



yo ... qp-' J 





1 



... 10 



(51) 
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To solve equation ()14|) . we leave Z an arbitrary j9 by p complex matrix, and plug into 
equation (|T^ our known form for X and Y: 



ZX - qXZ = - ^ hiY'-^ 

iy^O mod n 



(52) 



X 



Zu{l -q) ... Zij{q^ ^ -q) ... 
Ziiil - (f) . . . Zij{q^-^ -q') ... 



\ 



\ 



( 



...J 



\ 



hi ... 
&2 h ... 
h &2 '■• 



(53) 



Note that the sum truncates since = 0. These equations are solved by sets of z's such 
that 



hi 



Zij = for i + 1 < j 



(54) 
(55) 

(56) 



with the p — 1 elements -Zj^j+i still unrestricted. The other elements are uniquely fixed by 
the above recursion relation. 

To solve equation (fT5|l . we plug in our forms giving: 



YZ - qZY = - ttiX'-^ 

i^O mod n 



-qzi2 -qzn ■■■ 

2^11 — C[Z22 Zi2 — qZ23 . . . Zip 



\ 



^ Vl{x) 



Vliqx) 



V 



(57) 



(58) 



Analyzing (j3Hj) . we see that given the relationships between the Zij required to satisfy (dlj), 
all but the diagonal equations are automatically satisfied. Using the relation xV'{x) = 
f{x) — f{qx), with / defined as above, we get a system of p equations: 



-xqzi2 = -f{x) + f{qx) 
qx{zi2 - qz23) = -fiqx) + f{q^x) 

q^^'^x{zp^2,p-i - qzp-i,p) = -fiq^^'^x) + fiq^'^x) 
q'-'xzp^i,p = -fiq^^'x) + fiq^x) 



(59) 
(60) 
(61) 
(62) 
(63) 
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We have p equations with p — 1 variables, so the equations must not be hnearly independent 
to be solvable. Adding the equations together collapses nicely, giving: 

fix) = fiq^x) (64) 

Therefore only for values of x satisfying this equation will there exist a fractional brane 
solution with dimension p < n. 

We could have also taken Y diagonal in the first step, interchanging the role of X and Y. 
An identical derivation to the one above gives that solutions exist only at points such that 

g{y) = g{q^y), where y is the clement yn of Y. Thus wc have two different sets of fractional 
brane solutions for each p, of number approximately equal to the degree of the polynomial 
f{x) — f{q^x) plus the degree of g{y) — g{q^y)- For each of these either X"' or Y'^ are zero, 
this is or V vanish at that point. 

D. Singularities of the hypersurface equation 

We derived the hypersurface equation for bulk brane solutions. This was 

n—l n—1 

uvw = Hit - fiq'x)) + Hit - g{q'y)) - T (65) 

i=l i=l 

The singularities of this hypersurface equation are easy to locate, we search for points where 
all the partial derivatives vanish. The partial derivative with respect to w shows that either 
V or M is 0. Assuming v = and the requirement that the original hypersurface equation is 
solved fixes t to be equal to one of the f{q^x). Therefore the partial derivative with respect 
to t is satisfied only if f{x) — f{q^x) for some i (i.e. there is a double root in Q{t,u)). 
Likewise, if = we get singularities at points where g{y) — g{q^y). Note that unlike 
previous work where U{Y) = 0, we do not have any remaining codimension two lines of 
singularities, these have been deformed leaving behind only codimension three singularities. 

We now see that the fractional brane solutions directly correspond to the singularities 
found above. We derived that fractional brane solutions exist at points where f{x) — f{q^x) 
or g{y) — g{q^y) for some i, which exactly correspond to the singularities of the hypersurface 
equation. However, we have to be more careful. These are singularities of the Calabi-Yau 
geometry only if x or y are different from zero. This is because we have to take partial 
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derivatives with respect to u, v and not with respect to x, y to see if the point of the Calabi- 
yau is singular or not. Away from the branching locus of u = and v = ?/" calculating by 
taking partial derivatives with respect to x or y is allowed. At the origin of the branching 
locus one can actually calculate directly in the hypersurface equation and find that so long as 
ai or bi are non- vanishing, then this is a fake singularity. Also in that case one can not solve 
for the fractional brane solutions as we did above, as both X and Y will be nilpotent, and 
neither can be argued to be diagonalizable. Notice also that for each singularity there are two 
fractional branes associated with it. If f{x) = f{q^x) for i < n, then f{{q^x)) = f{q^~^{q^x)), 
and both of these solutions are associated to a brane sitting at the same value of u. A similar 
argument holds for v. Also, one can check that for generic values of the superpotential 
coefficients, all of these singularities lead to conifolds (double points). 



1. Relevant Quantum Deformations of the Hypersurface 

We know that when considering the full quantum system these singularities could be 
resolved by quantum deformations of the hypersurface equation . These are related 

to partial gaugino condensates ^. We will show that the number of moduli of these de- 
formations is exactly equal to the number of fractional brane solutions divided by two (i.e. 
the number of singularities). The idea is that for each of these conifolds singularities we 
have one relevant quantum deformation (we will also be considering marginally relevant) as 
defined in j^: Given a Calabi-Yau manifold and its holomorphic 3-form f2, a deformation 
of the form bif{u, v, w, t) is relevant iff 

lim 4- / ^ ~ C>(logA) (66) 

A^oo dbi Jb. 

where Bj are the set of non-compact cycles and A is a cutoff required by the non-compactness 
of the Bj. We can analyze what forms of deformation satisfy this by assigning charges and 
looking at the scaling behavior of the holomorphic 3-form. 

First, we assign charges to the variables so that {F{u,v,t,w) = 0), the hypersurface 
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equation, has charge 1. This fixes the degrees of the variables to be 

= - (67) 
n 

ru = \ (68) 

r. = - (69) 
m 

= 1 - 7 - - (70) 

where I is the degree of f{x) and m is the degree of giy). The holomorphic 3-form is 

^ du Adv A dw , 

^ = 9F (71) 

at 

and has weight Vu + + rw + rt — 1 — ^- If we have a deformation^ of the form F — biuH^, 
then the condition to be relevant is that 



d du A dv A dw 



dbi dF+F 

' at 



(72) 



bi=0 

scales as less then (this includes the marginal case) . Taylor expanding and taking the only 
surviving term we get that 

OF du Adv A dw 
'dt (^j^ 

which scales as 

+ + +_ + !_ 21-- 74 

/ n I m I m \ nj 

Requiring that this be < gives 

an + bl < l(n - 1) (75) 

The number of quantum moduli is therefore equal to the number of solutions (a, b) to 
this equation for given integers n,l,m.. We also have quantum deformations of the form 
F — bivH^, which are counted by the solutions to 

an + bm < m{n — 1) (76) 



^ This is the most general form since terms with both u and v can be absorbed into w. 
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2. Counting Singularities 

Here we will show that the number of moduli is equal to the number of singularities or 
fractional brane solutions (actually pairs of fractional branes). We found that there were 
singularities at points where f{x) = f{q^x) or gijj) = g{q^y). Since / is a degree / polynomial 
in X and g is degree m, f{x) — f{q^x) is also a degree / polynomial and has / — 1 solutions, 
except in the case where g*' = 1, where it has only / — 2 solutions (we aren't counting as a 
solution, we'll add that special case later). Also, we note that there is a symmetry between 
the solutions, since if there is a root xq that satisfies /(x) = f{q^x) then x'q = g*xo satisfies 
/(x) = /(g"~*x). Since Xq = x'q^ these points are the same in the full geometry associated 
to the M, V variables. Therefore we count the number of roots by 

1 — 2 li = mod n 
d{t) = { (77) 

/ — 1 otherwise 



and 



m — 2 mi = mod n 
e{i) = { (78) 

m — 1 otherwise 



as m 



The total number of singularities is then 



(79) 



2 



plus 1 for the singularity at 0, which is actually made from n identical fractional branes at 
that point. Now the counting question has been reduced to a purely mathematical question: 
Is g equal to the number of solutions of an + 6/ < l{n — 1) plus the number of solutions 
of an + bm < m{n — 1) (plus one to count 0)? The answer is yes, and although you can't 
perform the sum explicitly you can show that the sum counts the number of solutions of 
the inequalities exactly. Therefore, for each singularity there exists a pair of fractional 
branes that behave like the Klebanov-Strassler setup 0], and a gaugino condensate on these 
fractional branes resolve the associated singularity. This way we count one quantum moduli 
per conifold of the deformed Calabi-Yau. 
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E. Explicit Solutions for small n 

1. n = 2 Case 

Restricting to this case for illustrative purposes, we see that q — —1 and the above 
solutions become 

X^x\ ^ ° I (80) 

F = y I ° (81) 

Z-i (82) 

The center elements are: 

u^X'^^x'' (83) 

v = Y'^ = y'^ (84) 

w = -\ir{Z'^) ^ det Z ^ z'^ - zl- {z[f (85) 

t = xyz (86) 

(In this case, — |trZ^ = detZ, since Z is traceless.) Therefore 

uvw = x^y'^iz^ -zl- {z[f) (87) 

J.2 2 2 2 2 2/ /\2 foo\ 

-X y z^-x y [z^) (88) 



\j=l mod n / \j=l mod n / 

UVW = t'-'^(j2 (^i+2jU^ " ^ b^+2jv'^ (90) 

This is the hypersurface equation. We can see this is equal to the general form derived 
earlier by noting: 

uvw = t^ + f{x)f{-x)+g{y)g{-y) (91) 
^{t- f{x)){t - fi-x)) + {t- g{y)){t - g{-y)) - f (92) 
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since f{x) + /(— x) = in this case. 

Where are the singularities of this equation? First, the partial derivative with respect to 
t gives that t = 0, and the derivative with respect to w gives that uv = 0, therefore either 
u or V must be 0. The partial derivative with respect to u gives: 



vw 



(93) 



The partial derivative with respect to v gives similarly: 



uw 



(94) 



In the f = case, the first equation is solved by the roots of the rhs polynomial, and the 
second simply fixes w to a constant. If m = it's the other way around. So our singularities 
are at these points: 



t = 0, f = 0, Mj = Roots 



t = 0, M = 0, f j 



Roots 



vj=0 
0=0 



Wi 



4:Ui 



4vi 



(95) 



(96) 



Of course, looking closer at these polynomials, we see that (with generic values for the 
constants) the roots are exactly the roots of f{x) and g{y) squared as defined above. The 
other roots of (jH!^ and and not singularities since that solution does not satisfy the 
original hypersurface equation. 

The fraction brane solutions are trivial since they are one-dimensional. Solving the F- 
term equations: 

xy = (97) 

implies either x or ?/ is 0. Let y = 0, then: 

(98) 
(99) 

and 



2zx = —bl 

z = -^ 
^ 2x 



(100) 



i=0 



This is exactly the singular points where f = as derived from the algebraic geometry 
above. The a; = case is identical. 
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2. n = 3 Case 



The general solution to the equations is: 



Y = y 



1 

g 

V ^2 y 

'^O l\ 

1 
1 



^2 



z\ + z^q + zq^ 
z[q^ 
z'^q^ 



^2? 

Zi + zq + Z2(f' 
z[q 



\ 



The elements of the center are u — X^ — x^, v — Y^ — y^, and 



(101) 



(102) 



(103) 



w 



= -tr(Z^) = det Z = + + ^1 + {z[Y + (4)^ - ?,z{ziZ2 + ^i4) 



(104) 



(Again, in this case the trace is equal to the determinant expression, however that is no 
longer true in n > 4). As usual, t — xyz (we always arrange for it this way to simplify the 
hypersurface equation). The hypersurface equation is: 



uvw ^t^ + x^y^izf + zl + {z[f + {z'^f - 3z{ziZ2 + z'^z'^) 



-3'»IEi^)IEf?f )-3'MEi^)(i:f!? 



q I \ — i- — q 



(105) 
(106) 

(107) 



q I \ i — q^ 

j=o ^ J \j=o ^ 



(108) 



This actually factors as noted in the general solution, and can be written in a simpler form 
as: 

3 3 

UVW = Hit - f{xq')) + Hit - g{yq')) - (109) 
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(there are no cross terms in this case either), with / and g given in equations ()37p and ()38p. 

Let's derive this using the determinant form instead. We can write the matrix t — f[X) — 
g{qY) out exphcitly in terms of t and the functions / and g: 



det 



^ t — f{x) —xyz'2 —xyz[ ^ 

-xyz[ t - f{qx) -xyz'2 
y -xyz'2 -xyz'^ t-f\q^x) j 



(110) 



Writing out the terms we get: 
2 

- /(^'^)) - - - (3t - f{x) - f{qx) - f{q\))xVz'A (HI) 

1=0 

First we note that f{x) + f{qx) + f{q^x) = 0. It is also easy to prove that uv{z'iY + uv{z'2)^ = 
g{.y)g{.qy)9{q^y) from g{y) = xyz'^ + xyz'2, g{qy) = qxyz'^ + q^xyz'2, and g{q^y) = q^xyz'^ + 
qxyz'2. Therefore it simphfies to 

2 

- /(^'^)) - 9{y)9{qy)9{q'y) - stxVz'A. (112) 

We can also show directly by using the above that g{y)g{qy) + g{,qy)g{,q^y) + 9{q^y)9{y) = 

—?)x'^y'^z'^z'2., so we've derived explicitly that: 

2 2 
uvw = \[{t - fiq'x)) + Hit - g{q'y)) - (113) 

i=0 i=0 

as claimed. So as in the U{Y) = cases, there are singular points where f{x) = f{xq^) for 
p < n, and in addition we also have singular points at g{y) = g{yq^) as well. 

III. MATRIX MODEL 

We now use the correlation between this theory and a bosonic matrix model, following 
the procedure used in Il5|, and the combinatorial solutions to multi- matrix models worked in 

n 

|4fl| . Studying this related matrix model should allow us to solve for the quantum corrected 
moduli space. Our calculation above found the moduli space for a probe brane in the 
geometry, and we will split the eigenvalues of the matrices in the matrix model into sets that 
correspond to bulk branes, X type fractional branes, and Y type fraction brane eigenvalues. 
We use this gauge choice to integrate out cross term matrices and derive an effective potential 
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for the branes in the bulk. This is our quantum moduh space for probe branes as derived 
from the matrix model, and the derivation parallels exactly the one in [l5|. 
The associated matrix model has potential 



W{X, Y, Z) = Nfi-^ [tT{XYZ - qXZY) + tr V{X) + tr U{Y)] 
with n correlated eigenvalues corresponding to a single brane in the bulk. 

A. Quantum Deformations of Potentials 

To compute the quantum corrected moduli space, we follow the procedure 
divide the matrices as follows: 

\ / v; "v^ . n . . . n \ 



(114) 



in 0. 



We 



X 



X, 
Ai 



X 



bfy 











An. 



Xf^., 



X 



Y 



fy J 



••• 






a; 



v 



a: 



Zhb Zf^ Zf^fy 

y z fyb Zfyf^ z fy j 



ri y j 



;ii5) 



with Xf, as the n by n matrix corresponding to the brane in the bulk, eigenvalues for 
fractional branes on the X type singularities, and Uy eigenvalues for fractional brane on Y 
type singularities. The partial diagonalization of X and F is a gauge choice, and therefore 
gives you a Vandermonde determinant which gives a contribution to the effective potential 
of the form 

2 log(X6 - Ai) + log(Ai - A,) + 2 trfe(log(n - AD + X] loglA^ - A;.) (116) 

when you raise it up into the exponential. 

When you substitute this into the potential, it is easy to show that you get twice the 
number of terms as in and that Xbfy,Xfyb, Ybj^, Yj^b and all the off diagonal Z's appear 
quadratically and can be integrated out, giving more contributions to the effective potential 



5^tr,(log(X,-gA,))-5^tr,(log(X,-g-iA,))-5^tr,(log(n-gA:))-5^tr,(log(n-g-'A:)) 

i i it 

(117) 
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where we are ommiting terms that don't include the bulk brane matrices. So in total the 
effective potential for the bulk brane in this model is: 



W = iV/i-i(tr((X,n^fe - qX,Z,Y,) + V{X^) + U{Y^) 



N ^ 

i 

N ^ 



2 log(X5 - A,) - log(X, - gA,) - ^ log(X, - q-'X 

i 

2 log(n - K) - ^og{Y, - q\[) - J2 log(n - q-'K 



(118) 
(119) 

(120) 



At this point we simply have two copies of the results of |la], and can apply the results 
of that paper. The quantum version of f{x) is the same as in that case: 

n-l 

f{X) = f{X) + J](l - q^)X--^R,{u) (121) 



where 



^ oo 



(A"^+^') 



k=0 ' 



fc+1 



(122) 



are the generating functions of the moments of the eigenvalues of X. Likewise we now also 
have the quantum version of g(Y): 

n-l 

g{Y) = g{Y) + ^^(1 - q^)Y^^'m'^{v) (123) 
i=i 



where R' is defined with the A'. 



B. Quantum Relations 

The point of using the matrix model to derive the quantum moduli space is that we claim 
these are equal to the deformations allowed by holomorphy we derived for the hypersurface 
equation. In fact we derived that the hypersurface equation is modified to: 

Q{t, u, v) = r + e-\0,{u) + 0[{v)) + e-^{02{u) + O'^iv)) + ■■■ (124) 
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where the O's and 0"s are polynomials in u and v respectively^. From the definitions given 
above, we can derive relations like: 

Oi{u) = J2f{qW{q^x) (125) 

i<j 

-02{u) = fiqWiq'^mq'^) (126) 

i<j<k 

and so on for the rest of the moduli. These are the quantum relations for the matrix model, 
and there are also the corresponding equations for 0'^{v) as well. Since these are identical 



to the quantum relations in jl5|, we can use the results that these equations give recursion 
relations for the moments of the eigenvalues (A'^). For example, 

^a^(A'+"^) = ^ (127) 

i 

comes from requiring that the coefficients of the non-polynomial terms of the RHS of the 
first quantum relation vanish, as they must since Oi{u) is polynomial in u. The second 
and higher quantum relations give similar recursion relations between the moments. The 
identical equations also hold with hi and A ^ A'. 



C. Loop Equations 

The quantum relations that we derived above from relating the matrix model to the 
q deformed field field theory should also be encoded directly in the loop equations of the 
matrix model. The loop equations have the universal form of 

Y.^^^f^)^f-) = {^^<i^) (128) 

where the sum is over all monomials /i,/2 such that /i0*/2 = /(0)- Here we are using a 
simplified notation where: 

(/(0)) = (^tr(/(0))^ (129) 

^ The terms proportional to can be eliminated by linear shifts in t, or can also be seen to vanish from 
the scaling argument given above. 
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1. First Equation 

By taking the variations 5Z = X^'ZX^-'^, SZ = Y"'ZY''-"', SX = X''+^ and SY = Y''+\ 
the loop equations we get are 

= (YX'^ZX'"'^+^) - q{YX"'+^ZX''~"') (130) 

= (xy""+izr'=-™) - g(xr'"zr'=-™+^) (131) 

k 

fi J2 (^"') {X''-'^) - {V\X)X''+^) = {X^+^YZ) - q{X^+^ZY) (132) 

m=0 
k 

^ ^^ym^^yfc-m^ - {W {Y)Y''+'^) = {Y'^+'^ZX) - q{Y^+^XZ) (133) 
m=0 

It seems we have k + 2 equations for the k + 2 variables (yx^^ZX^"""^^^) and another + 2 
equations for the k + 2 variables l^XY"^ ZY^^"^^^) which determines them in terms of the 
{X^) and (y™-) respectively. However when k = Q mod n, these equations are no longer 
linearly independent, so for the equations to even have a solution, we sum the first k + 1 
equations weighted by q: 

k 

= ^g'"(rX'"ZX*^-"+^) - g(rX'"+^ZX^-'")) (134) 



m=0 



- q^-'{YX''+'Z) (135) 



k 



fi J2 {X"") (^'"'") - {V\X)X''+^) (136) 



m=0 

^m\ I vk—m\ 



Since {X"^){X is symmetric over the sum from to k, we get that: 

(K'(X)X^+^) = ^ ^ /i(2 - g*" - g-") (X"^) (X'^"™) (137) 
which is exactly the quantum relation p27p derived above. Clearly equations ()13ip and 



()133p will give the v quantum relation. All the comments given in 
untwisted one-point functions still apply. 



about twisted versus 



2. Higher Equations 

The second and higher quantum relations should be obtained in the matrix model loop 
equations from higher term variations like 6Z = X"^ Z'^Y X^^^~^ . This is essentially the 
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same technique as used in jl^, however we no longer can set words with unequal numbers 
of y's and Z's to zero, since we no longer have the U{1) charge symmetry between them. 
We will look at variations and use the same idea as above. 
First, we take the variation 6Z = X"'Z^YX''+^-"': 

fi{X"'Z){YX^+'^'"') + = {X"'Z^YX''+^-'^Y)~q{X"'+^Z^YX^+^-'^Y) 

(138) 

This is again k + 2 equations for A; + 3 unknowns (the lower terms have already been solved 
for in the first order equations in terms of the (X"^)), so we need another equation to solve 
the system. Actually it will be convieniant to have two extra equations (with one more 
unknown), which come from the variations 6X = ZYX^^"^ and 5X = YX^^'^Z: 

k+l 

/i ^(ZFX™)(X^+i-") - (ZrX^'+V(X)) = (rZZFX'=+2) - q{ZY ZY X''+'^) (139) 

m=0 
fe+1 

11 (FX'") {X^+^-'^Z) - {YX^+'^ZV'{X)) = {ZYZYX^+^) - q{ZZYYX''+^) (140) 

m=0 

Again in the case where k = mod n, these equations are linearly independent and we get 
a consistency condition by summing them. By summing p38|l you can rewrite in terms of 
the sum of p39|l and p4()|l . and you get the constraint: 

k+l k+1 

- g'")/i(X™)(ZrX'=+^-'") + J](g - g^'")^(X^'"Z)(FX'=+i-™) = 

m=0 m=0 

(ZFX'=+V(X)) + g(FX^+'Z\/'(X)) (141) 
But we can see from the variation 6Z = X^ that (yX^) = 0, so the equation simplifes to: 

k+l 

^(1 - g'")/i(X™)(ZyX^+^-'") = (ZrX^+V(X)) + q{YX^+'^ZV'{X)) (142) 

Now this equation has reduced to essentially the same equation as in [13]. You can solve the 
{ZYX^) terms from the first order equations in terms of the (X'^), and that determines a 
recursion relation between them. The modification of the above to derive the second order 
Y quantum relation is clear. 
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